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Abstract 

t> 

We study the fluxes of anti-deuterons that could be produced by anni- 
hilations in the galactic halo of Dark Matter particles with multi-TeV 
mass and a large annihilation cross section, as indicated by the re- 
cent PAMELA results. The model of Minimal Dark Matter (MDM) 
is an example in this category. We find that the fluxes are well within 
1 the reach of planned experiments for DM candidates that annihilate 

mainly into quark pairs, and also extend into the multi-GeV range 
^ above the background. They are instead suppressed and concentrated 

at too large energies to emerge above the background if the main an- 
nihilation channel is into gauge bosons (such as in particular W + W~ 



O 

in the MDM case). 

> 

b 1 Introduction 

03 

Cosmology and astrophysics provide several convincing evidences of the existence of Dark Mat- 
ter [1]. The observation that some mass is missing to explain the dynamics of galactic or cluster 
rotations dates back to the '30s [2]. The observations from weak lensing [3], for instance in the 
spectacular case of the so-called 'bullet cluster' [4], provide evidence that there is mass where 
nothing is optically seen. More generally, global fits to a number of cosmological datasets 
(Cosmic Microwave Background, Large Scale Structure and also Type la Supernovae) allow to 
determine very precisely the amount of DM in the global energy-matter content of the Universe 
at tt DM h 2 = 0.110 ±0.005 [5]. 

All these signals pertain to the gravitational effects of Dark Matter of the cosmological and 
extragalactical scale. Searches for explicit manifestation of the DM particles that are supposed 



1 



to constitute the halo of our own galaxy have instead so far been giving negative results, but 
this might be on the point of changing. 

Indeed, a recent flurry of results from a number of cosmic ray experiments has shown a few 
surprising anomalies that can be interpreted as the first indirect signatures of annihilations of 
DM particles in the galactic halo, on top of the standard cosmic rays (CR) of astrophysical 
origin. 1 We review below the precise indications of these experiments, but we anticipate that 
the central piece of our motivations are the PAMELA results, that indicate the possibility 
of a heavy DM particle (multi-TeV) that can annihilate with a large cross section (of the 
order of 10 _23 cm 3 /sec or more, much larger than the expected thermal cross section of ~ 
3 • 10~ 26 cm 3 /sec) into quarks or SM gauge bosons. 

In indirect DM searches, the key point is to look for channels and ranges of energy where it is 
possible to beat the background from ordinary astrophysical processes. This is for instance the 
basic reason why searches focus on fluxes of anti-particles, much less abundant in the Universe 
than the corresponding particles. The proposal to look in particular for fluxes of anti-deuterons 
in the cosmic rays was put forward for the first time in ref. [10]. The basic observation was 
that anti-deuterons are rarely produced in standard processes. Especially the low energy range 
(below about 1 GeV) is kinematically disfavored for astrophysical production and therefore 
open to the possible emergence of an exotic contribution like the one from DM. More recently, 
anti-deuterons from DM annihilations have been studied in [11] and [12]. 

Motivated, as said, by the PAMELA results, we here focus our analysis on the case of DM 
particles that are heavy (multi-TeV masses) and have a large galactic annihilation cross section. 
We compute and present the fluxes of anti-deuterons adopting a model independent approach, 
i.e. considering different possible masses and primary annihilation channels. But we also study 
the specific predictions of the model of Minimal Dark Matter [29] , that naturally falls into this 
class. 

The rest of the paper is organized as follows. In sec. 2 we describe the DM models on 
which we concentrate our attention and their motivation on the basis of the data. In sec. 3 
we describe the production of anti-deuteron from DM annihilations and its propagation in the 
galaxy, discussing in particular the astrophysical parameters and the associated uncertainties. 
In sec. 4 we present our results, comparing them with the expected sensitivities of planned 
experiments. In sec. 5 we draw the conclusions. 

2 The Dark Matter framework 

The recent results from the PAMELA satellite [13] show 

• a significant excess above the expected smooth astrophysical background and a steep rise 
of the positron fraction e + /(e + + e~) above 10 GeV up to 100 GeV [14], compatibly with 
previous less certain hints from HEAT [15] and AMS-01 [16]; 

• no excess in the p/p energy spectrum [17] compared with the predicted background, up 
to the maximal probed energy of about 100 GeV. 

If interpreted in terms of Dark Matter annihilations, these results indicate [18] either 

x It is entirely possible that these anomalies will be explained in terms of a better understanding of standard 
CR [6] or rather peculiar astrophysics, such as one or more local pulsars [7], nearby sources of CR in galactic 
spiral arms [8] or exploding stars [9] . We are of course instead interested in the possibility that DM is at the 
origin of the signals. 
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Table 1: Masses and annihilation cross sections of a heavy DM particle that fits the PAMELA 
positron excess (see [18]), rather independently from the annihilation channel. The 1 TeV case is 
quite disfavored by the antiproton results and is kept mainly for reference purposes. The last line 
refers the Minimal Dark Matter model reviewed in sec. 2.1. 

(i) a DM particle x °f an y mass (above about 100 GeV) that annihilates only into leptons 
(xx ~ > e + e~, fi + fi ~, r + r ~), not producing therefore unseen antiprotons or 

(ii) a DM particle with a mass around or above a few TeV, that can annihilate into any channel 
(i.e. XX W + W~ , ZZ, bb, ti, light quark pairs and the leptonic channels above) possibly 
producing anti-proton fluxes at energies above those currently probed by PAMELA. 

We will of course concentrate on the second possibility in this paper, as annihilations into 
leptonic channels do not produce any anti-nucleons and a fortiori no anti-deuterium. Instead, 
annihilations into gauge bosons or quarks copiously produce anti-proton and anti-neutron fluxes 
and therefore entail a possibly large production of anti-deuterons. We will consider Dark Matter 
masses Mdm equal and above 1 TeV as indicated in table 1, being aware that the smallest value 
(1 TeV) is kept mainly for reference purposes, as it is already excluded by the non-observation 
of anti-protons. 

The PAMELA results also require a very large annihilation cross section, in order to match the 
large observed flux in positrons. Based on the analysis of [18] we adopt the benchmark values 
indicated in table 1. These very large cross sections, much larger than the typical annihilation 
cross section required by DM thermal production in cosmology, ~ 3 ■ 10~ 26 cm 3 /sec [1], can be 
justified in specific models in terms of some enhancement mechanism which is effective today but 
not in the early universe (such as a resonance [18, 19] or Sommerfeld [20, 30, 21] enhancement, 
the presence of an astrophysical boost factor due to DM substructures -unlikely [22]-, or a 
combination of these). For the purposes of our model independent approach, they are taken as 
an input required by data. 

In addition to the PAMELA data, the balloon-borne experiments ATIC [23] and PPB- 
BETS [24] have recently reported results indicating the presence of an anomalous peak in the 
spectrum of cosmic ray e + + e~ at about 500-800 GeV. The HESS telescope has also published 
data [25] in the range of energy from 600 GeV up to a few TeV, showing a steepening of the 
spectrum which is compatible both with the ATIC peak (which cannot however be fully tested) 
and with a feature power law with index —3.05 ± 0.02 and a cutoff at ~2 TeV. If interpreted in 
terms of DM annihilations, the ATIC peak would clearly pin down the DM mass at about 1 TeV, 
and as a consequence select the possibility (i) above, making the generation of anti-deuteron 
fluxes impossible. However, the interpretations of the results from these experiments have to 
rely on MonteCarlo simulations that are tested only up to much smaller energies. Moreover, 
the agreement between the different balloon experiments is not optimal, e.g. the EC data [26] 
fail to clearly indicate the presence of the peak. We work therefore under the assumption that 
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the possibility (ii) above is the correct one. The upcoming data from the FERMI satellite [27] 
will be crucial to settle the issue. 

If the signals in PAMELA are due to DM annihilations, one should also worry about the 
associated emission of high energy gamma rays from the galactic center (GC), the galactic ridge 
and the Sagittarius dwarf galaxy, and also of radio waves from the galactic center (produced 
by synchrotron radiation in the strong magnetic field by the electrons and positrons from 
DM annihilations). These can impose stringent constraints on the annihilation cross section. A 
complete, model independent calculation has been carried out in this respect in [28]. It is found 
that, for most annihilation channels, the cross sections required by PAMELA are excluded if 
the DM distribution in the galaxies follows a benchmark NFW profile (reviewed below). But 
choosing a smoother profile and/or assuming that a part of the cross section is due to an 
astrophysical boost factor that would not be present in dwarf galaxies and the Galactic Center 
due to tidal disruption re- allows the cross sections of table 1. How realistic these choices and 
assumptions are will be clarified by further observations, studies (especially of dwarf galaxies 
DM distributions) and simulations. 

In summary: the recent PAMELA data select a heavy, multi-TeV DM particle with a large 
annihilation cross section (see table 1) as a possible candidate for explaining the excess in 
positrons and the null result in antiproton fluxes (the tension with the ballon experiment data 
in the e + + e~ spectra has to be settled by further data and we do not consider it in the 
following, while the bounds from 7 rays and radio signals can be avoided with assumptions on 
the DM distribution that will be tested in the future). 



2.1 Minimal Dark Matter 

A particular model featuring DM particles of the kind favoured by the PAMELA data is the 
Minimal Dark Matter (MDM) model [29, 30, 31, 32, 33]. It consists of a minimalistic approach 
to solve the dark matter problem, not inspired by more ambitious Beyond-the-Standard- Model 
ideas such as supersymmetry or extra dimensions. It is constructed by adding to the SM one 
new n-plet of the SU(2) L gauge group and assigning quantum numbers such that a viable dark 
matter candidate is found: a stable neutral particle still allowed by current DM searches. To 
automatically insure stability, n-plets that possess decay modes into SM particles consistent 
with renormalizability are excluded. To keep it really minimalistic, scalar n-plets with quartic 
couplings to the SM higgs are rejected as well. The only interactions with SM particles are then 
of gauge type and the only new parameter is the tree level mass M MDM of the new multiplet. 
This mass is fixed by the requirement that the thermal relic abundance equals the measured DM 
abundance. Direct searches for DM exclude all n-tuplets with Y 7^ as they have interactions 
with the Z boson that lead to an elastic scattering cross section with nucleons which is 2 - 
3 orders of magnitude above present bounds from direct DM searches. All in all, a fermionic 
quintuplet with Y = and mass M MDM = 9.6 TeV is singled out as the most successful 
candidate. 

MDM particles annihilate into W + W~ at tree level and into ZZ, 7Z and 77 at one-loop. Due 
to the Sommerfeld enhancement, the cross sections for these annihilations are much higher than 
what is normally expected for WIMPs. One finds [30]: 

( av ) ww = 1.1 • HT 23 — , (av) zz = 3.2 • 10~ 24 — , (<rv) yZ = 6.5 ■ 10~ 25 — . (1) 

sec sec sec 

Note that this is an example of a model in which the annihilation cross section is of the typical 
weak order of 3-10 _26 cm 3 /sec at thermal freeze out, and thus produces the right DM abundance 
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Qdm — 6 ■ 10~ 27 ^- (c - ^) -1 (the so-called "WIMP miracle"), but then acquires the larger values 
of eq. (1) at later times. This is because the Sommerfeld enhancement is roughly inversely 
proportional to the particle velocity [30] and the particles in the early universe had much 
higher velocities (/? ~ 0.1) than in the galactic halo today (/3 ~ 10~ 3 ). 
If a boost factor due to DM overdensities in the halo 

B ~ 50 (2) 

is assumed 2 , the predictions of the MDM model for an excess of positrons in cosmic rays due to 
DM annihilations are well matched by the PAMELA results (see addendum to [31] and [33]). 



3 Anti-deuteron production and propagation in the galaxy 

3.1 Primary fluxes from DM annihilations 

We consider the annihilation of a dark matter particle x with its antiparticle x m t° pahs of 
bosons and quarks: 

XX - W + W~, jZ, ZZ, cc, bb, ti, (3) 

where XX ~ * °c stands representative for annihilations into light quarks. Of course, other 
annihilation channels, such as XX ~ y 77 are a l so conceivable but would not produce any p and 
n and thus no d. 

Extending the SM to incorporate DM particles, we simulate these processes using the Mad- 
Graph 4.4 Monte Carlo event generator [34] for the masses of the DM particles x indicated in 
table 1. We then hand over the resulting bosons and quarks to PYTHIA 8.1 [35] for showering 
and hadronization, explicitely requesting the n not to decay, as they are needed for the forma- 
tion of d. The produced p-fluxes as a function of the kinetic energy T agree very well with the 
fluxes presented in [36] 3 . The fluxes of ra-fluxes are essentially equal to the p ones, as expected 
from isospin invariance. 

3.2 Anti-deuteron formation by coalescence 

Anti-deuterons are then produced via the fusion process of a p-n pair. This is described by the 
very simple coalescence model [37]. In this model the so-called "coalescence function" ^(kp, k n ) 
describes the probability for a p-n pair to yield by fusion an antideuteron. It actually depends 
only on the difference of the momenta kp — k n = 2 A and is strongly peaked around A = 0: 

kp^kn^^- (4) 

2 We stick to these given values of av and B for definiteness in our analysis of MDM, but we note that such 
values can vary by about one order of magnitude -thus possibly reducing B to about 5 and correspondingly 
increasing av- within the 3 a range of J^dm (see [31, 33]). A larger av exacerbates the tension with gamma 
ray constraints from the galactic center and dwarf galaxies, but a smaller boost would be preferable for the 
agreement with numerical simulations. 

3 Our fluxes are of course a factor of 2 lower as we keep n from decaying. We note that the fluxes in their 
fig. 3 most probably refer to dN/dx and not dN/dT as indicated. 
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x=r/M DM x=r/M DM 

Figure 1: Spectra of antideuterons per annihilation from DM annihilations in the indicated primary 
channels, for the masses of table 1. 



The ^-density in momentum space is thus written as the p-density (in momentum space) times 
the probability to find an n within a sphere of radius po around kp. 

d 3 Nj 4tt , d 3 N n d 3 N, 

J-density in probability to find p-density in 

momentum space «■ within a sphere ^ momentum space 
of radius p around k p 
in momentum space 

where 7 is the Lorentz factor and the "coalescence momentum" p is a free parameter con- 
strained by data on hadronic production [58, 59, 60]. The central value preferred in [51, 12] is 
Po = 79 MeV, while a window of 66 - 105 MeV is given in [12]. It is energy independent for a 
wide range of energies [51, 58]. 

In brief, one assumes that the antinucleons merge whenever they are within a momentum po of 
each other [56, 57, 12]. After some algebra, one gets the following formula for the <i-spectrum 
at the site of production, where T denotes the kinetic energy: 

dT d 3 Po m n m p Jf2 ~ 2m _ T _ V dT p j Tf= S \ dT n ) Tn= T _A 

Notice that we apply the coalescence prescription to the fluxes of p and n as produced by 
the MonteCarlo code. A more precise computation, that we do not perform, would require to 
coalesce the p and n pair event per event within the MonteCarlo itself. 

Fig. 1 shows our results for the primary antideuteron fluxes. In fact, all quark channels lead 
to fluxes that are quite similar among each other, as do respectively the gauge boson channels. 
We therefore select here and in the following as exemplary the cases of the bb and the W + W~ 
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channels. One sees that the gauge boson channels lead to fluxes that are less abundant and 
concentrated at high energies. This allows to foresee a less promising signal (after propagation) 
with respect to the one from primary quarks. 



3.3 Propagation 

Once the antideuterons have been produced, in each given point of the DM halo, they propagate 
in the galactic region and eventually reach the Earth. In this section we review the processes 
to which they undergo during their trip and the formalism that allows to compute the resulting 
final spectra. 

Assuming steady state, the diffusion-loss equation describing the change in the differential 
number density J\f(r, z,T) = 4^ of cosmic ray particles when travelling through the galaxy in 
cylindrical coordinates r, z is [44, 45, 46, 47, 48] 

V[K(T)VN{r,z,T)-V c {r,z)N{r,z,T) ] - T(T)H{r, z, T) + £{r,z,T) 



(spatial) diffusion convection annihilation source 

d 



| V-V c {r ,z) d 



dT 



fV^T) 



<)T [(b ion (T)+b Co »i(T))N(r,z,T)] 



adiabatic loss due to galactic wind ionization and Coulomb losses 

d 



(7) 



dT 

s 



1 -±^K pp (T)M(r,z,T) 



d_ 
+ df 



(3 2 K pp (T)^Af(r,z,T) 



first & second order reacceleration 

The last four terms are energy redistribution terms. They can be neglected [38] for the fluxes of 
primary anti-deuterons. The basic physics reason is that all these processes happen in the thin 
disk of the galaxy, where the anti-deuterons spend little time (they rarely cross the galactic 
plane). We keep however the annihilation terms that have the more drastic effect of removing 
a particle from the flux 4 . The same approximation is used in [36], where the procedure to 
re-include those terms is also briefly discussed. 

To be able to analytically solve the diffusion equation (7) one has to employ a simplified 
description of the spatial distribution of matter in the Galaxy. The so-called two-zone disk-halo 
model is adopted, where the Galaxy is described as a thin gaseous disk (of radius R = 20 kpc 
and half-height h = 100 pc) embedded in a thick, cylindrical diffusive halo. The half-height 
L of the halo is a free parameter of the model, constrained from the study of Boron/Carbon 
(B/C) CR ratios [44, 49]. 

Spatial diffusion arises because charged particles interact with the galactic magnetic field 
inhomogenities. Diffusion is an energy-dependent process as higher energy particles probe larger 
spatial scales than lower energy ones. The energy-dependent diffusion coefficient is given by 
[44, 48] 

K{E) = K f3n 5 (8) 

where (3 as usual denotes the velocity in units of c and the rigidity, 1Z = ^ for a particle of 
atomic number Z. Kq and S are again constrained by B/C data. 

The presence of a galactic wind V c of cosmic rays pointing outwards from the galactic plane is a 



4 This result applies to p as well. Notice that positrons, also used for the indirect search for Dark Matter, 
have energy loss terms corresponding to Inverse Compton and bremsstrahlung: these effects are not relevant 
for d (and p) since 3> m e + . 
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Figure 2: Cross sections for the different processes of interaction of light nuclei on the interstellar 
gas, and the experimental datasets on which they are based. 



well established fact. This leads to convection and an energy loss due to the adiabatic expansion 
of the plasma (neglected here), as described in the diffusion equation (7). The strength of this 
wind, V c , is again a parameter to be determined by fitting to B/C-data. 

The transport parameters L, K , S, V c providing the maximal, median and minimum primary 
antideuteron flux compatible with B/C analysis are reported in table 2. We will later study 
the effect on the propagated fluxes of choosing one or the other of the different sets. 

The third term in eq. (7) describes the annihilation of the anti-deuterons with H and He in 
the interstellar medium (ISM) in the thin galactic disk. T(T) is therefore defined as: 



where 



T(T) = 2M(^)r ann (T) 
\T) = (n H + 4§n H , 



V °meA T ), 



(9) 

(10) 



with nu « 1/cm 3 the hydrogen density and nn e ~ 0.07 nn the Helium density. In order to 
calculate 
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corresponds to anti-deuterons interacting inelastically with the ISM but sur- 
viving the collision, loosing only a fraction of their energy, one needs the total, elastic and 
non-annihilating inelastic cross sections for the process d + H [12]. For what concerns the total 
cross-section, no experimental data exist, but data for the charge conjugate reaction dp can be 
obtained from the PDG [50]. We assume therefore 



non— ann 
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dp 
tot 
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dp 
tof 



(12) 



As for the elastic cross section, no data are available even for the charge conjugate reaction. 
We therefore make the assumption 

<f « 2o?, (13) 



and data on a^f can again be obtained from the PDG. At energies > 10 GeV, we make this 
assumption for the total cross section as well. As it would be complicated to calculate the 
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cross section cr™" -8,1111 = a (dp — > dX) due to the complex interplay of the underlying NN cross 
sections and selection rules, it is obtained from experimental values (of the symmetric system 
pd — > Xd as there are - again - no experimental values for the reaction dp — > dX). o~(pd — > Xd) 
is found by summing up the pd — > (nir)pd cross sections experimentally available. Channels 
not known experimentally are assumed to give only negligibly small contributions [51]. Figure 
2 illustrates the different ingredients that we have just described. 

The source term contains a primary term (the signal from DM), a secondary term (the 
astrophysical background) and a tertiary term 

&{r, z, T) = Q prim (r, z, T) + Q sec (r, z, T) + Q tert (r, z, T). (14) 

The tertiary source term merely describes the fact that some antideuterons move to higher 
energy bins by interaction with the ISM, others to lower energy bins. As already stated above 
we neglect this energy redistribution. We will take the propagated background fluxes from the 
literature (see below), so we only solve the equation for the primary source 



The sum is over annihilation channels with p, n in the final state and annihilation cross section 
of the DM particles (av) k and dN^/dT are the source spectra computed in sec. 3.1. The 
coefficient rj depends on the DM particle being or not self-conjugate: for example, for fermions 
it is 1/2 for Majorana and 1/4 for Dirac particles. For definiteness (and also for consistency 
with the specific case of Minimal Dark Matter) we choose to focus on the Majorana case. The 
Dirac case just amounts to rescaling our results by a factor 2. 

For the distribution of DM in the galaxy, p DM (r,z), we consider, as customary, three different 
models determined by numerical simulations. The Navarro-Frenck-White profile [52] 

Pnfw(0 = Ps— f 1 + — ) , (16) 



the Einasto profile [53, 54] 

PEinasto(r) = ps exp 

and the isothermal profile [55] 
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The values for the parameters r s and p s of the three models are given in table 3. While they 
sensibly differ at the Galactic Center (from the cored isothermal profile to the more cuspy 
NFW), they do not at the proximity of the Earth. We will later study the effect on the 
propagated fluxes of changing the choice of profile. 

Neglecting energy gain and loss terms, using cylindrical coordinates and Bessel-expanding 
Af(r, z, T) and Q prim (r, z, T) 

oo 

Af(r,z,T) = J2 N ^ T )Jo (19) 

i=l 



9 



case 


L [kpc] 
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kpc 2 
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V c [*f] 


DM halo model 


r s in kpc 


Ps in ^ 


max 


15 


0.46 


0.0765 
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NFW [52] 


20 


0.26 


med 
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0.7 


0.0112 


12 


Einasto [53, 54] 


20 


0.06 


min 
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0.85 


0.0016 


13.5 


Isothermal [55] 


5 


1.16 



Table 2: Transport parameters providing the Table 3: Parameters of the density profiles 
maximal, median and minimal primary an- for the Milky Way DM halo, 
tideuteron flux compatible with B/C data [12] 



where is the ith zero of Jo(x), the equation to solve in the halo is [44, 45, 46, 47, 48] 

d 2 V d £ 2 \ O prim (z T) 



dz 2 K(T) 3z B?J iy ' J K(T) 

Imposing the boundary condition Ni(z = L) =0 and ensuring continuity with the solution in 
the disk, where the equation to solve is 

d 2 V d f 2 \ 2hT ann (T) 

^ * N t (z, T) = ' j ^)^(0> T) (21) 



dz 2 K(T)dz R 2 ) iV ' ' K(T) 

(there are no sources in the thin disk, but there is annihilation on the ISM, described by 
r ann (T)), one gets the solution 5 
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Finally, for the flux $>(r,z,T) = j-J\f(r,z,T) near the Earth one therefore gets 



$(r = r Sun ,2; = 0,T) = — cWl 



B L m\ 



4tt y (T + m d -) 2 
X b^[rmf) L ) A sinh (\SiL) \ 



(26) 



i=l 

where -B is the boost factor. 

The last step needed is to include the average solar modulation effect: the solar wind 
decreases the kinetic energy T of charged cosmic rays such that the flux $Earth (denoted as 
TOA, for Top Of the Atmosphere) of anti-deuterons that reach the Earth with kinetic energy 
^Earth and momentum pEarth is approximately related to their flux in the interstellar medium 
as [61] 

<£Ea rth = ^f^, T = T Earth + |Ze|0 Fisk , p 2 = 2m,T + T 2 . (27) 



5 With 

Vc. , , \ . . fS 



y t (z)=2 I e*v[^{z- z')Uinhl^{z- z'))QY™{z')dz' (22) 



A l = 2hT*™ + V c + KS l -coth[^) (23) S . = J^_ +A ^ (24 ) 
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M DM = 10. TeV, DM profile: Einasto, <cr v>=7xl0 -22 cm 3 s"' M DM = 10. TeV, propagation: med, <cr v>=7xl0 -22 cm 3 s"' 




J in GeV/n T in GeV/n 

Figure 3: Left panel; Propagated (top- of-the- atmosphere (TOA)) d flux, including solar modulation, 
from a DM with mass 10 TeV annihilating into bb, fixing the DM halo profile at Einasto and changing 
the propagation parameters (min, med, max). Right panel: Fixing the propagation parameters at med 
and changing the halo profile (NFW, Isothermal, Einasto,). 



In this effective formalism the so-called Fisk potential 0Fisk parameterizes the kinetic energy 
loss. We assume a value of </>Fisk = 0.5 GV, characteristic of a minimum of the solar cyclic 
activity, corresponding to the period in which most of the observations were done in the second 
half of the 90's and are being done now. 

4 Results and discussion 

Infig.s 3- 5 we show our results. Fig.s 3 and 4 illustrate the impact of the different astrophys- 
ical choices for one sample DM of Mdm = 10 TeV annihilating into bb, while fig. 5 summarizes 
our main results for all the considered DM candidates and also reports the indication of the 
current experimental results and future sensitivities. 

The background estimation reported in all our results is taken from the most recent calcu- 
lation of [12], fig. 1, which agrees and slightly updates the detailed analysis in [51]. Background 
d are produced by the spallation of high-energy cosmic ray (CR) p, He and p on the interstellar 
(IS) gas made of H and He in the galactic disk: the dominant source is CKp on IS H for energies 
above 1 GeV, CR p on IS H and He for lower energies. Tertiary antideuteron production is 
important in this case. Notice the characteristic shape that decreases at E < few GeV, opening 
an interesting window for the emergence of possible DM signals, as anticipated in the Intro- 
duction. For definiteness, the propagation scheme of the secondary d is fixed at med and the 
solar modulation at solar minimum (although, in principle, one should adapt the propagation 
parameters to the ones considered for the primaries). The uncertainty on such background is 
quoted at about a factor of 2 in [51]. 
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TOA rf-flux from XX ~* bb TOA d— flux from XX ~* ^ 

M DM = 10. TeV, Einasto, med, <tr v)=7xl0 -22 cmV M DM = 10. TeV, Einasto, med, (a v>=7xl0" 22 cmV 




T in GeV/n T in GeV/n 

Figure 4: Left panel; Propagated d flux from a DM with mass 10 TeV annihilating into bb, fixing 
an Einasto DM profile and med propagation parameters, with and without the solar modulation ef- 
fect. Right panel: The same flux (including solar modulation) for different choices of the coalescence 
momentum pq. 



As it can be inferred from fig. 3a, the choice of propagation parameters has quite an impor- 
tant impact on the predicted fluxes: the differences can amount to almost 2 orders of magnitude 
at the lowest energies. This is consistent with the analogous results in the case of anti-proton 
fluxes (see e.g. [31]). 

The choice of halo profile, instead, has a very limited effect, as shown in fig. 3b. This is un- 
derstood in terms of the fact that essentially no antideuterons come from as far as the region 
of the Galactic Center where the profiles differ significantly. This result is also in analogy with 
the antiproton case. 

The solar modulation effect, shown in fig. 4a, has a certain effect in the lowest energy 
portion of the spectrum (it amounts to about a factor of 2 below 1 GeV), while the choice of 
the coalescence momentum within the intervall quoted in sec. 3.2 has of course a simple overall 
normalization effect that amounts to about a factor of 4. 

All in all, it is evident that the uncertainties ascribed to astrophysics and nuclear physics 
are very important, spanning a few orders of magnitude in the normalization of the fluxes (the 
shapes are instead more robust). Future improvements on B/C CR measurements might allow 
to restrict the ranges for the propagation parameters, a dominant source of uncertainty for the 
predictions. 

Searches for anti-deuteron fluxes have been performed by BESS [39], a balloon-borne mag- 
netic spectrometer experiment that collected data between 1997 and 2000. The null result 
imposes an upper limit on the flux in the explored energy range at 

0f ESS < 1.9 • 10" 4 (m 2 sec sr GeV/n)" 1 , 0.17 GeV/n < T~ d < 1.15 GeV/n (95% C.L.). 

(28) 
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TO A d-{\ux from XX ~^ bb or W + W (Einasto & med propagation parameters) 



current BESS limit 




T in GeV/n 



Figure 5: The final fluxes of anti-deuterons for the different heavy DM candidates discussed in the 
paper and for Minimal Dark Matter compared to the background (gray shaded area) and to expected 
experiment sensitivities. Here we have assumed an Einasto DM profile, med propagation parameters 
and included the solar modulation effect. 



Future experimental capabilities rely on the performances of the AMS-02 cosmic ray experiment, 
to be installed on the International Space Station, and the dedicated GAPS balloon-borne 
experiment. AMS-02 [41, 42] is expected to be sensitive to a flux 

tf — > 4.5 • 10" (m 2 sec Sr G.V/n)- { % <j$- < g < ; 8 GeV/n (2g) 

GAPS [40] will use a dedicated technique: an incoming anti-deuteron is slowed down in the 
detector and forms an exotic atom with one of the nuclei of the detector (replacing a shell 
electron); eventually the atom de-excites emitting several hard X-rays and the captured an- 
tideuteron annihilates on the nucleus producing a characteristic hadronic shower correlated with 
the X-ray signature. The Long Duration Balloon (LDB), Ultra Long Duration Balloon (ULDB) 
and possible satellite-borne (SB) flights foresee a sensitivity respectively to a flux [40, 11] 

^GAPS-LBD > L g . 10 _ 7 (m 2 sec gr Q eV / n yl 
^GAPS-ULBD > 3 q . 10 - 8 ( ffi 2 gec gr Q e y /n) -l 

^gaps-sb > 2 6 . 10 - 9 ^ m 2 gec gr GeV/n)- 1 o.l GeV/n < T s < 0.4 GeV/n (31) 

The final fluxes of antideuterons at Earth, after solar modulation and for the benchmark 
choice of astrophysical parameters (Einasto, med) are presented in fig. 5 for the masses, annihi- 
lation channels and corresponding annihilation cross sections (table 1) that we discussed. The 



0.1 GeV/n < T s < 0.25 GeV/n (30) 
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Figure 6: 'Best' (DM profile: NFW, propagation parameters: max) and 'worst' ('isothermal, mm) 
scenarios for the detection of heavy dark matter via anti-deuterons. We do not include the 1 TeV case 
here as it is already excluded by the non- observation of anti-protons. We do not include fluxes for 
annihilation to W + W~ either as they are below the background anyway. 



fluxes from primary W + W~ are a few orders of magnitude below the background and also con- 
centrated at high energy, so essentially undetectable. These features can be essentially traced 
back to the characteristics of the primary fluxes presented in fig. 1 . The case of Minimal Dark 
Matter (green line in fig. 5) falls in this particular category. On the contrary, the fluxes from 
primary bb are large (and relatively flat) at low energies, where the sensitivity of the upcoming 
experiments is located. Fluxes from larger mass candidates are slightly more suppressed (the 
source term depends on M^, only partly compensated by a larger cross section). It is inter- 
esting to note that some of the fluxes can extend into the multi-GeV region (possibly probed 
by the AMS-02 mission) while remaining above the background. 

It is worth reminding that different choices of the astrophysical parameters lead to order of 
magnitude changes in the overall normalization, as discussed above: for instance for max prop- 
agation parameters and NFW profile the predicted fluxes (from bb) for any DM mass lie well 
within the AMS-02 area ('best case scenario'), while for min and isothermal only the 5 TeV 
and 10 TeV candidates are predicted to produce a flux that skims through the GAPS region 
('worst case scenario'). This is illustrated in fig. 6. 



5 Conclusions 

We have studied the fluxes of antideuterons produced by heavy Dark Matter annihilations in 
the galactic halo. We have proceeded in a model independent way, considering a range of DM 
masses and corresponding large annihilation cross sections motivated by the PAMELA results 
on the excess in positrons, if interpreted in terms of Dark Matter. These are shown in table 1. 



14 



We have also considered the specific model of Minimal Dark Matter. 

We have computed the primary fluxes of p and n at production via the PYTHIA MonteCarlo 
code, generated the d fluxes (presented in fig. 1) via the coalescence analytical model and then 
propagated them across the galactic halo with the semi-analytic method of the two-zone model. 

Fig. 3 and fig. 4 show the impact of the different choices for the astrophysical propagation 
parameters, halo DM profile, the coalescence momentum parameter and of the solar modulation 
effect. The uncertainty on the final predicted fluxes can be large, even of one or two orders of 
magnitude. 

Fig. 5 shows our main results: the fluxes of antideuterons at Earth for the different masses 
and primary annihilation channels, together with the predicted experimental sensitivity. For 
the case of annihilation into a bb pair (exemplar of any annihilation into a quark pair) the fluxes 
are large, concentrated at low energies and therefore in the optimal conditions for detection 
in the upcoming experiments. Fluxes from the W + W~ channel (gauge bosons in general) are 
instead not relevant, and this in particular applies to the model of Minimal Dark Matter that 
has W + W~ as the dominant channel. 

Our results can be used for predicting the antideuteron fluxes in any other given precise 
model, where branching ratios for the different primary annihilation channels are known, simply 
linearly combining and rescaling appropriately the spectra that we presented. 

In summary, we find that if the PAMELA signals in positrons are correctly pointing to 
DM particles that are heavy (multi-TeV) and with very large annihilation cross sections, then 
the associated fluxes of antideuterons from the same DM anihilations are very promising for 
the upcoming experiments (GAPS, AMS-02), if the dominant annihilation channels are into 
quarks and not gauge bosons (and, of course, not pure leptons). Astrophysical uncertainties 
have however a sizable impact. 
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